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Abstract. We study the quantum operation of coupled superconducting flux qubits 
under a microwave irradiation. The flux qubits can be described as magnetic 
dipole moments in the limit of weak microwave field amplitude consistent with usual 
experimental situations. With the Hamiltonian for coupled qubits under a microwave 
field, we show that a strong coupling enables to realize the high performance controUed- 
NOT gate operation. For practical quantum computing we analyze the effect of 
microwave on switching function of phase-coupled qubits. 
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1. Introduction 

Superconducting Joseplison junction qubits are one of the most promising candidates for 
implementing quantum computation because macroscopic coherence of superconductor 
is robust against noises from environment. Recent experiments for superconducting 
charge [1] and flux [21 [3] qubits have reported much longer coherence times. For practical 
quantum computing also the high performance coupled-qubit operations need to be 
achieved. For flux qubits an XY-type of coupling between two qubits has been achieved, 
where only a SWAP gate operation has been demonstrated [Ij. The controUed-NOT 
(CNOT) gate, as the basic element of the universal gate provides the simplest design 
for scalable quantum computing. By using an Ising-type of coupling between two flux 
qubits, a CNOT gate operation has been experimentally demonstrated for inductively 
coupled flux qubits [6] with rather weak coupling. 

In this study we obtain the Hamiltonian for coupled three- Josephson-junction 
qubits under a microwave field in weak microwave amplitude limit, consistent with 
usual experimental situations. The Hamiltonian is written in terms of the magnetic 
dipole moment of qubit which constitutes the qubit-microwave coupling constant in 
experiment f7] and phenomenological Hamiltonians [8]. For a proper parameter regime 
for the CNOT gate operation the Hamiltonian is reduced to a block-diagonalized form. 
Each diagonal part of the Hamiltonian corresponds to different control qubit state. The 
discriminating oscillations of occupation probability of coupled qubits give rise to the 
CNOT gate operation. With the obtained Hamiltonian for coupled qubits we show 
that strongly coupled qubits can achieve the high performance CNOT gate operation. 
The fluctuation effects of the microwave field as well as the static magnetic field are 
discussed. We consider, for example, phase-coupled flux qubits to achieve a strong 
coupling. For practical quantum computing the phase-coupling scheme of flux qubits 
provides a switching function. We discuss the effects of microwave on the switching 
operation of phase-coupled qubits. 

2. Hamiltonian of coupled flux qubits interacting with a microwave field 

When only a static flux is penetrating a flux qubit, the fluxoid quantization condition 
for the qubit loop is given by 2?™ + 27r/ — 01 — 02 — 03 = with integer n. The 
reduced flux is denoted as / = ^st/^o with the static external flux $st threading 
the qubit loop and the superconducting unit flux quantum $o = h/2e. Using the 
relation 0i = 27r(n + /) — 02 — 03, the energy levels of the qubit is written as 
Es{f) = Eji [1 - cos(27r(?2 + /) - 02,s - 03,s)] + Ej2{l - cos 02, J + ^j3(l - COS 03,,), whcre 
02(3), s with s G {i, t} is the value of 02(3) for the state |s) at the local minima of Ues 
and depends on /. Here | |) (| ])) is the diamagnetic (paramagnetic) current state of 
the flux qubit. 

For the Rabi oscillation of a flux qubit, a microwave field with frequency u is 
applied; /<^(t) = ($mw/'^'o) coscut, where $mw = BS with the microwave magnetic field 
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B and the area of the qubit loop S. Using the fluxoid quantization condition 

27rn + 27r(/ + /^(t))-0i-02-03 = O, (1) 

the effective potential of the qubit is written as Ues{4>) = Eji[l — cos(27r(/ + f^jit)) — 
02 — ^s)] +-Ej2(l — cos + -Ej3(l — cos ^s) with n = 0. Normally, in usual experiments, 
the applied microwave fields are in the range of <l>mw ^ *^'o- Then, the effective potential 
is written as 

t/eff(0) ~^ji[l-cos(27r/-02-03)] (2) 

+ Ej2{l - COS 02) + Ej3{l - COS 03) 

+ 27rEji($mw/$o) costutsin(27r/ - 02 - 03), 

where 02(3), s can be assumed to be constant in time due to •I'mw ^ '^'o- 

In the basis of the qubit current sates {|t) , one can obtain the qubit 
Hamiltonian. The diagonal components E^^^^ of the qubit Hamiltonian consist of the 
static Eg and oscillating parts; E^,s{t, f) — Es{f) — MgB cos ut with 

M, = -Sh sin(27r/ - 02,. - 03,s) = SI, (3) 

where the qubit loop has the circulating current / = — /csin0i with /c = 2ttEji/^q and 
the magnetic dipole moment of the qubit is Mg. Hence, in terms of the magnetic moment 
of the qubit loop and the interaction with the microwave field, the qubit Hamiltonian 
is rewritten as 

n = Ho-M-B{t), (4) 

where Ho = E^ |i) (i| |T) (T| -t^di) (T| + IT) (i|) and tg is the tunneling amplitude 
between the two states in the qubit, which comes from the charging energy of the 
Josephson junctions [9j. The interaction between the magnetic dipole moment and the 
microwave field is described by 

M ■ B(t) = MB cosa;t(|i) (i| - |T) (T|), (5) 

where we set M| ^ — M| = M. Here note that, although the microwave field just 
threads the qubit loop, not applied on the qubit directly, the dipole magnetic moment 
of qubit and the microwave interact with each other through the fluxoid quantization 
condition of Eq. ([T]) . 

For two coupled flux qubits we consider that the left qubit is the control qubit. In 
this case the flux J'l is adjusted far away from the degeneracy point so that the tunneling 
process t^ in the left qubit is negligible, i.e., ti/lEis — Ei^s\ ^0. As a consequence, the 
two-qubit Hamiltonian becomes block-diagonalized. Hence the problem is reduced to 
that of two independent qubits under a microwave irradiation, 

H=J2[EsAfL, fn) - {Mls + MR,,)Bcosujt]\s, s'){s, s'\ 

s,s' 

-tR\s,s'){s,-s'\, (6) 

where Ess'^JljIr) is the energy level of coupled qubits, Ml(H)s is the magnetic dipole 
moment of the left (right) qubit in |s) state, and —s denotes the opposite pseudo-spin 
state of s G {i, T}- Here we set Mr^ ^ -Mr-^ = Mr and Mli ^ -Ml^ = Ml- 
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To clearly describe the Rabi-type oscillations, we employ rotated coordinates for 
the coupled-qubit states as |00) = cos(^;/2)| H) + sm{ej2)\ it), |01) = -sin(^|/2)| U 
) + cos(ei/2)| it), |10) = cos(^t/2)| Ti) + sin(^t/2)| TT), and |11) = -sin(^t/2)| Ti 
) + cos(6'|/2)| tl) with tan6's = 2tfi/\Esi — Esi\. Let us consider, for example, the case 
that the control qubit states is | J,). Then we see the relations 

lii) (ill - lit) (itl = cos^i(|00)(00| - |01)(01|) (7) 

-sin^i(|01)(00| + |00)(01|), 
lii) (111 + lit) an =|00)(00| + |01)(01|. (8) 

If 2tR/\Esi — Egi^l = 0, i.e., sin 6*1 = 0, the off-diagonal terms in Eq. ([7]) does not appear 
so that the transition between the qubit states does not occur even for a resonant 
microwave field. Hence the tunneling tji between the states | ii) and | it) plays a 
key role in responding to the microwave field. For the control qubit state | |) we also 
perform the transformation and obtain the two-qubit Hamiltonian, 

^ = E [^po{t)\pO){pO\+£p^{t)\pl){pl\ 

p=0,l 

+ a,MnBcostut{\pO){pl\ + |pl)(pO|)], (9) 

where 

^PP'W=C'(/L,/ii) - [{-iYMl + i-iy'pMBcosut 

(10) 

with ao = sin6'|, ai = sin 6^1, Pq = cos9i, and Pi = cos9p 

The Hamiltonian in Eq. iQ is valid in the weak microwave amplitude limit 
'^'mw ^ *&o- However, to perform a gate operation via a Rabi oscillation in experiments, 
^jaw should satisfy a more strict condition that the Rabi frequency fl is much smaller 
than the energy gap uq, i.e., 

n = MB/h<^uo. (11) 

In this regime, the rotating wave approximation (RWA) can be applied and a well- 
behaved Rabi-type oscillation can be observed. In other words, the applied microwave 
field should be in the range of = BS ^ Hujq/I. From the experimental parameters 
for the flux qubits in Ref. [10], this condition reads <I>mw ^ 10^^<I>o- For the Rabi 
frequency Q/27i ^ 600 MHz which coincides with usual experimental situations, we 
find that the corresponding amplitude ~ 10~^$o provide a well-behaved Rabi-type 
oscillation in this study. 



3. Controlled-NOT gate operation 

In this section, we consider a concrete system, for example, two phase-coupled qubits 
[m [121 [HI HH [15] in Fig. [H From the Hamiltonian of Eq. we show that the 
CNOT gate operation can be achieved with a high performance for a strong coupling. 
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In a previous study [TT] the CNOT gate operation was analyzed without microwave 
irradiation. The energy levels Sp^/^fi, of coupled qubits are shown in Figs. [2](a) and 
I C 1 cLS cL function of fn with fixed fi = 0.49. 

At degeneracy points (solid lines) in Figs. [2](a) and (c), = = Eq, E-^^ < E-^-^, 
9 1 = tt/2 and 9^ = tan~^(2ti^/|£^|| — E-^il). At this point we have the diagonal elements, 
Sqo = —tji — MiB cos ujt + Eq, £qi = tji — MlB cos u!t + Eq, and the constant energy gap 

uq = £Qi{t) - SQoit) = 2tR. (12) 

On the other hand, the off-diagonal term aQM^B cosut (p = 0) with cto = 1 gives rise 
to a dynamical evolution between the states |00) and |01), i.e., a Rabi-type oscillation. 
In Figs. [2] (b) and (d) the occupation probabilities Ppp' of \pp') states during the Rabi- 
type oscillations are shown as a function of time when the initial state is prepared as 
V'(0) = (|00) + |10))/V2. 

As shown in Fig. [2], the microwave field with the resonance frequency u = uq 
induces the Rabi oscillation between the states |00) and |01) owing to the off-diagonal 
term with p = in Eq. ([9]). For the weak coupling case of Fig. Mh) the states |10) 
and 1 11) also oscillate simultaneously in response to the microwave field through the off- 
diagonal term aiMi^S cosa;i:(|10)(ll| + |11)(10|) with p = 1, while they are stationary 
for the strong coupling case in Fig. EJ^d). 

For a discriminating Rabi oscillation, the coupling strength J should be larger than 
the tunneling rate, J > t^. The coupling strength is given as J = {E^i^ — E^i)/A at 
the degeneracy point where E^^ = E^^- [HI |12j. Further, as discussed in Eq. ffTTj) . for 
Rabi-type oscillations the energy gap ujq = 2tR should be much larger than the Rabi 
frequency Q = M^B/h. Consequently, for a high performance of CNOT gate operation, 
we see the criteria 

n < c^o(= 2tR) < J. (13) 
In Fig. [2](a), at the degeneracy point fn ~ 0.4994, the tunneling amplitude is greater 




Figure 1. Phase-coupled two flux qubits connected by a loop which has two dc- 
SQUIDs providing a switchable coupling between flux qubits. The Josephson coupling 
energy of a junction in the dc-SQUIDs with phase difference 4)a(b)i is E'j. The dc- 
SQUIDs have threading flux fa{b) which can be adjusted to control the coupling 
between two flux qubits. When fa = fb = 0, 3. dc-SQUID can be simplified as a 
single Josephson junction with the Josephson coupling energy 2Ej. 
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than the couphng strength, i.e., tji/h ^ 2GHz and J/h ^ 0.6GHz. Thus, the oscillations 
are not discriminative. To improve the discrimination of oscillations, one need to increase 
J larger than the value of tR. 

From the energy levels in Fig. [2]^c), increasing J is shown to make the distance 
farther between the degeneracy point at ?a 0.4945 and the other degeneracy point. 
This property makes it possible to perform a discriminative Rabi-type oscillation. In 
Fig. El^d), we plot the oscillations of occupation probabilities Ppp/ as a function of time 
by using the Hamiltonian in Eq. ([9]). Figure EJ^d) shows that the microwave field at the 
resonance frequency u = ujq induces the Rabi oscillation between the states |00) and 
|01) driven by the off-diagonal term for p = in Eq. ([9]), while the states |10) and |11) 
(p = 1) do not respond to the microwave field. From Eqs. ffTOj) and (|T2|) we have the 
relation 



001= Suit) -£io{t) = \J{AJf+ujl + 2(3iMRB cosut, (14) 

which shows that the difference \u!i — cuqI between the energy gaps increases as the 
coupling strength J increases at the degeneracy point. Here the the latter time- 
dependent term is negligibly small compared with J due to the criteria of Eq. f|T3|) . 
Hence the states |10) and |11) hardly respond to the microwave with frequency ujq, 
resulting in the CNOT gate operation via discriminative Rabi-type oscillations at 
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Figure 2. Energy levels E^p, of coupled qubits in Fig. [T] for (a) a weak coupling 
J /h =0.6GHz and (c) a strong coupling J /h = 5GHz and occupation probabilities 
of coupled flux qubits during Rabi-type oscillations at the degeneracy point (b) 
fa « 0.4994 and (d) Jr ~ 0.4945 for the weak and strong coupling case, respectively. 
Here we choose the parameters as = 0.49, Eji/Ej — 0.75, and tij//i=2GHz. The 
initial state is ipiO) = (|00) + |10))/V2 and the Rabi frequency is n/2'K = 600MHz. 
Ess' with s,s' e {i,T} are shown as thin dotted hues in (a) and (c). At Qt = (odd)7r 
the CNOT gate operation is expected to be achieved. 
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nt={odd)7i in Fig. (SJ^d). 

Now we discuss the effect of ffuctuations of microwave field as well as the static one 
on the CNOT gate operation. During the two-qubit operation both the static flux fL{R) 
and microwave flux fu}(t) give rise to the noises, Sf^R^it) and Sf^(t), which destroy the 
qubit coherence. These noises appear in the fluxoid quantization condition, 

-(t>L{R)l - 4>L{R)2 - 4>L{R)3 = 0. (15) 

Let us first discuss the noises in the left qubit. Combining the fluctuation of static flux 
6 flit) into that of the microwave field Sfi^{t), the net effect of both noises is expressed 
as a random fluctuation 5Bi{t) in the microwave field fuj{t) = {BS/^q) cosut. As a 
result, the magnetic field of microwave in the diagonal terms of Eq. (Q is rewritten as 
Ml{B cos ut + SBiit)). Since the diagonal term Ml{B cos cut + SB L(t)) does not change 
the energy gaps Uq and ui in Eqs. (fT2l) and ([HI), the fluctuation 5Bi{t) has no effect on 
the operations. Hence the fluctuations of both the static and microwave fields threading 
the left qubit loop are negligible. 

On the other hand, the fluctuations of fluxes threading the right qubit may cause 
the decoherence in the qubit states. Let us discuss the terms with p = and p = 1 in Eq. 
([9]) separately. For the p = terms in Eq. ([9]) the fluctuation of microwave flux, 6f^(t), 
can be combined into the fluctuation of static flux, SfR{t), in the fluxoid quantization 
condition for right qubit in Eq. f|T5l) . Since the degeneracy point is optimally biased 
with respect to the static flux f^, the first order fluctuation effect of both SfR(t) and 
^fuit) on qubit state dephasing will vanish at this point. 

For the terms with p = 1 the net fluctuation can be given by Mr{B cos ujt + 6Bji{t)) 
in the diagonal and off-diagonal terms. Since 9-^ = tan~^(tR/2| J|) ^ and thus ai ~ 
for a sufficiently strong coupling, the off-diagonal term with p = 1 will not appear. Thus 
the fluctuation SBnit) in the off-diagonal terms hardly gives rise to dissipation of qubit 
states for a sufficiently strong coupling, but the fluctuations in the diagonal terms may 
cause dephasing. 

4. Effect of a microwave on the switching function of coupled qubits 

For practical quantum computing, the quantum operations should be manipulated by 
a switchable coupling. To discuss this, in the model of Fig. [1], we introduce two dc- 
SQUIDs in the connecting loop of the phase coupled qubits. Here the external fluxes 
fa and fb vary from zero to 0.5. When the fluxes threading the dc-SQUID loops are 
set as /a = /f, = (switch on), a dc-SQUID can be simplified as a single Josephson 
junction with the Josephson coupling energy 2Ej. The design of Fig. [1] is similar to 
that in a previous study [12]. The difference is just the direction of pseudo-spin of the 
right qubit. The present design is symmetric so that it may be more appropriate for 
real experiments. 

The current flowing in the connecting loop /' also gives rise to magnetic 
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moments in the dc-SQUID loops. /' depends on the two-qubit states such as /' = 
=F(27r/$o)2-EjCOS7r/'sin7r/^g = ±/q for the states, | |J.) and | |"f), respectively and 
otherwise I' = [12]. Here we set fa = fb = /'• Then the Josephson energy of the 
connecting loop U'jj ^{(f)ai,4>bi) = Yl'i=i[Ej{l — cos (j)ai)+E — cos (pbi)] has an additional 
oscillating term, 

2E'j27r{BS'/%) cos Trf sin vr/^g cos cut, (16) 

where we used the fluxoid quantization conditions for the dc-SQUID loops and the 
connecting loop, S' is the area of the dc-SQUIDs, and /^g = (0li + 0i?i)/27r. 
As a consequence, the connecting loop energy of coupled states becomes 

^L',. = ^L'-^L'5cosa;t, (17) 

where = -M', = M' and = = 0. M' = J^S" can be interpreted as the 
magnetic moment of the control loops. This magnetic moment arises by the interaction 
between the magnetic flux and the flowing current via the fluxoid quantization conditions 
of the dc-SQUID loops. If we include these terms in the Hamiltonian of Eq. (Q, the 
net effect is just the shift of magnetic moments, 

MsL ^ MsL + 0.5M', M,n ^ M^r + 0.5M', (18) 

remaining the physics qualitatively the same. 

The microwave field threading the dc-SQUID loops also generates noises. The 
noise from the microwave field can be introduced as a fluctuation of magnetic field 
0.5M'{B cosut + 6B'{t)) as before. 6B'{t) can be combined with the previous 
fluctuations 6B^ji){t) through the relation of Eq. (ITSl) . which does not generate 
qualitatively different effect. The static flux /' controls the coupling between two qubits; 
when /' = (0.5), the coupling is switched on (off). The static flux /' also generates 
noises in the dc-SQUID loops, but the switch-on (off) point, /' = (0.5), is an optimal 
point with respect to /' [12]. 



5. Discussions and summary 

For a direct comparison of CNOT gate operation performance as the coupling strength 
varies, the fidelity and concurrence [T6] oscillations are plotted as a function of time in 
Figs, and (b). The fidelity is calculated by the definition 

F{t) = Tr(p(t)pcNOT)/4, (19) 

where Pcnot is the matrix for the perfect CNOT operation and p{t) is the truth table 
amplitude at time t [6j. At t = 0, p(0) is the 4x4 identity matrix and the fidelity has 
a finite value F{0) = 0.5. The CNOT operation changes an initial product state into 
a maximally entangled state. Thus the maximal entanglement as well as the maximal 
fidelity corresponds to the perfect CNOT gate operation. Figure [3] shows that the 
deviations of the fidelity and concurrence oscillations diminish as the coupling strength 
increases. 
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Figure 3. (a) Fidelity and (b) concurrence oscillations with the initial state 
ip{0) = (|00) + |10))/\/2 during the CNOT gate operation for coupled flux qubits. 
The dashed (J = 0.6GHz) and solid (J = 5GHz) lines correspond to the weak and 
strong coupling of Fig. [21 respectively. At fit = Stt the coupling is switched off with 

fa^fb^ 0.5. 

At the end of two-qubit operations the phase-couphng is switched off with fa = 
ff, = 0.5. Since the magnetic dc pulse for switching-off has a finite rising time, the 
phase of qubit state evolves during the time, but this phase evolution is controllable 
by manipulating other parameters. The states of qubits can be detected by shifting 
the magnetic pulse adiabatically [3]. At the degeneracy point in Fig. [2|, the averaged 
current of qubit states vanishes. Thus, one can apply a finite dc magnetic pulse to shift 
the qubits slightly away from the degeneracy point to detect the qubit current states. 

In summary, we study the operation of two coupled flux qubits under a microwave 
irradiation. The flux qubits interact with the threading microwave field by the fiuxoid 
quantization condition and can be treated as magnetic moments for a weak magnetic 
field threading the qubit loop. By using the coupled qubit Hamiltonian we show that 
for a strong coupling the microwave-driven CNOT gate can be realized with a high 
fidelity. The strong coupling between flux qubits is obtained by a phase- coupling 
scheme. Introducing the switching SQUIDs in the phase-coupled qubit results in the 
renomalization of the magnetic moments of both qubits. The fluctuation effects of both 
the static and microwave fields are discussed. 
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